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Oh- Abstract 



We study the statistical properties of the triplet (a x ,a y ,a z ) of 
Pauli matrices going through a sequence of noisy channels, modeled by 
the repetition of a general, trace-preserving, completely positive map. 
We show a non-commutative central limit theorem for the distribution 
of this triplet, which shows up a 3-dimensional Brownian motion in 
qq ■ the limit with a non-trivial covariance matrix. We also prove a large 

deviation principle associated to this convergence, with an explicit 
\ rate function depending on the stationary state of the noisy channel. 

1 Introduction 

o 

In quantum information theory one of the most important question is to un- 
derstand and to control the way a quantum bit is modified when transmitted 
through a quantum channel. It is well-known that realistic transmission 
channels are not perfect and that they distort the quantum bit they trans- 
mit. This transformation of the quantum state is represented by the action 
of a completely positive map. These are the so-called noisy channels. 

The purpose of this article is to study the action of the repetition of 
a general completely positive map on basic observables. Physically, this 
model can be thought of as the sequence of transformations of small identical 
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pieces of noisy channels on a qubit. It can also be thought of as a discrete 
approximation of the more realistic model of a quantum bit going through a 
semigroup of completely positive maps (a Lindblad semigroup). 

As basic observables, we consider the triplet (a x , <7 y ,cr z ) of Pauli matrices. 
Under the repeated action of the completely positive map, they behave as a 
3-dimensional quantum random walk. The aim of this article is to study the 
statistical properties of this quantum random walk. 

Indeed, for any initial density matrix p in , we study the statistical prop- 
erties of the empirical average of the Pauli matrices in the successive states 
§ n (p in ),n > where $ is some completely positive and trace-preserving map 
describing our quantum channel. Quantum Bernoulli random walks studied 
by Biane in [TJ corresponds to the case where $ is the identity map. Biane pQ 
proved an invariance principle for this quantum random walk when pi n = \l- 

This article is organized as follows. In section two we describe the physi- 
cal and mathematical setup. In section three we establish a functional central 
limit theorem for the empirical average of the quantum random walk associ- 
ated to the Pauli matrices generalizing Biane's result pQ. This central limit 
theorem involves a 3-dimensional Brownian motion in the limit, whose co- 
variance matrix is non-trivial and depends explicitly on the stationary state 
of the noisy channel. In section four, we apply our central limit theorem to 
some explicit cases, in particular to the King- Ruskai-Szarek- Werner repre- 
sentation of completely positive and trace-preserving maps in M2(C). This 
allows us to compute the limit Brownian motion for the most well-know 
quantum channels: the depolarizing channel, the phase-damping channel, 
the amplitude-damping channel. Finally, in the last section, a large devia- 
tion principle for the empirical average is proved. 

2 Model and notations 

Let M 2 (C) be the set of 2 x 2 matrices with complex coefficients. The set of 
2x2 self-adjoint matrices forms a four dimensional real vector subspace of 
M 2 (C). A convenient basis B is given by the following matrices 

1 = ( I 1 ) a *= ( 1 I ) a v={°i ) <r ' = { 1 -1 ) 

where a x , a y , a z are the traditional Pauli matrices, they satisfy the commu- 
tation relations: [cr a .,cr y ] = 2ia z , and those obtained by cyclic permutations 



2 



of a x , <Ty, a z . A state on M 2 (C) is given by a density matrix (i.e. a positive 
semi-definite matrix with trace one) which we will suppose to be of the form 



a j3 
(3 1-a 

where < a < 1 and \(3\ 2 < a(l — a). The noise coming from interactions 
between the qubit states and the environment is represented by the action 
of a completely positive and trace-preserving map $ : M 2 (C) — > M 2 (C). 

Let Mi, A7 2 , . . . , .\ //,. ... be infinitely many copies of M 2 (C). For each 
given state p, we consider the algebra 

M p = M 1 <g> A7 2 ® . . . ® M k <g> . . . 

where the product is taken in the sense of W^*-algebra with respect to the 
product state 

uj = p® <g> $ 2 (p) ® . . . <g> $ fe (p) ® .... 

Our main hypothesis is the following. We assume that for any state p, the 
sequance $ n (p) converges to a stationary state p^, which we write as 



Poo 
12 



^oo 

^oo 1 -« 



oo 

'oo 



where < a:^ < 1 and |/3oo| < o^l — a^). 
Put 

= 2 Re(/3oo), v 2 = -2 Im(/3oo), w 3 = 2a 00 - 1. 
For every k > 1, we define 

x k = I (g) . . . ® 7 ® (a x — v i 7) (g) 7 (g) . . . 

y fe = 7 <g> . . . <g> 7 <g> (trj, - w 2 7) <g> 7 <g> . . . 

Z k = I®...®I®(<Tz-V 3 I)®I®... 

where each (cr — v, 7) appears on the k th place. 
For every n > 1, put 

n n n 

-^n = F„ = -^n = z k 

k=l k=l k=l 



with initial conditions 

-^o = Y = Zq = 0. 

The integer part of a real t is denoted by [t] . To each process we associate a 
continuous time normalized process denoted by 

= n-^X [nt] , Y} n) = n-^Y [nt] , = n^Z [nt] . 



3 A central limit theorem 



The aim of our article is to study the asymptotical properties of the quantum 
process (X^ , Y^ , z[ n ^ ) when n goes to infinity. This process being truly 
non-commutative, there is no hope to obtain an asymptotic behaviour in the 
classical sense. 

For any polynomial P = P(Xi, X2, . . . , X m ) of m variables, we denote by 
P the totally symmetrized polynomial of P obtained by symmetrizing each 
monomial in the following way: 

1 x ^ 

XnX i2 • • • Xi k > — Xi a(1) . . . x icr(k) 

where is the group of permutations of {1, . . . , k}. 
Theorem 3.1. Assume that 



(A) 



$> n {p)= Po o + o{—). 



Then, for any polynomial P of 3m variables, for any (t 1; . . . ,t m ) such that 
< ti < t 2 < . . . < t m , the following convergence holds: 



lim w 

= E 



?(!) d(2) ^(3) R (l) r(2) r(3) 



P(Bi;',B^,B^,....D-.D] : '.D] 
is a thi 

with covariance matrix Ct, where 



where (B^\ B^) t >o is a three-dimensional centered Brownian motion 



\-v\ -v x v 2 -v x v z 
C = I -viv 2 l—vl -v 2 v 3 

-V1V3 -V 2 V 3 l-vl 
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Remark : 

Theorem 13.11 has to be compared with the quantum central limit theorem 
obtained in [5] and [9]. In our case, the state under which the convergence 
holds does not need to be an infinite tensor product of states. We also 
give here a functional version of the central limit theorem. Finally, in j5] 
(see Remark 3 p. 131), the limit is described as a so-called quasi-free state in 
quantum mechanics. We prove in Theorem 13.1 1 that the limit is real Gaussian 
for the class of totally symmetrized polynomials. 
Proof: 

Let m > 1 and (to, ti, . . . , t m ) such that t = < ti < t 2 < . . . < 
t m . The polynomial P{X%\ Y™, Z%\ . . . , y£> , zf}) can be rewrit- 
ten as a polynomial function Q of the increments: xj:™\ Y^\ Z^\ X^' — 

y (n) v (n) _ v (n) y(n) y (n) y(n) _ y(n) y(n)_y(n) yW _ 7W 

A monomial of Q is a product of the form A^ . . . Ai k for some distinct 
ix, . . . , %k in {1, . . . , m} where A* is a product depending only on the incre- 
ments X t . — x[ n ^ , Y^ — Y^j l7 Zp — Zp_ . Since the Aj's are commuting 
variables, the totally symmetrized polynomial of the monomial A ix . . . A ik 
is equal to the product A ix . . . A ik . Consequently, it is enough to prove the 
theorem for any polynomial A*. 

Let i > 1 fixed, for every u±, u 2 , v 3 G R, we begin by determining the asymp- 
totic distribution of the linear combination 

[v\ + v\ + v\r^ [mxP - xgj + ^ 2 (Yf } - + «*(z£*> - (i) 

which can be rewritten as 

[ntj / 



Consider the matrix 
A -. 



1 






+ vl 




1 






+ z/| 





vi(a x - vj) + v 2 (a y - v 2 I) + v 3 (a z - v 3 I)j 

-v\V\ - v 2 v 2 + v 3 {\ - v 3 ) v x - iv 2 
vi + i v 2 -v\V\ - v 2 v 2 - z/ 3 (l + v 3 ) 



which we denote by 



a\ a 3 
a 3 a 2 
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ai,a 2 G R, a 3 G C. 

From assumption (A) we can write, for every n > 

{P) yp^ + faiS) l- aoo + n (4) 

where each sequence (4> n (i)) n satisfies: 4> n {i) = o{l/y/n). 

Let k > 1, the expectation and the variance of A in the state $ fc (p) are 
respectively equal to 

Trace(A$ fe (p)) 

and 

Trace(A 2 $ fc (p)) - Trace(A$ fe (p)) 2 . 
If both following conditions are satisfied: 

[nti] 

Trace(A$ fc (p)) = o(y/n) (2) 

k=[nU-i]+l 



and 



[ nt i, 

lim - V [Trace(A 2 $ fe (p)) - Trace (A$ fc (p)) 2 ] = a(t, -fc-i), (3) 

ra— >+oo 77, * — ' 

fe=[nti-i]+l 



then (see Theorem 2.8.42 in [3]) the asymptotic distribution of ([T]) is the 
Normal distribution A/"(0, a(tj — a > 0. 

Let us first prove (j2J). For every fc > 1, a simple computation gives 

Trace (A$ fc (p)) = [aia oo + a 3 /3 oo + a 3 (3 O0 + a 2 {l-a 0O ) + o(l/y/n)] = o(l/y/n), 

hence 

Trace(A$ fe (p)) = J] o(l/Vn) = o(y/n). 

k=[nti— i]+l fe=[ntj_i]+l 

This gives (J2]). 

Let us prove ([3]). Note that the sequence (Trace (A<I> n (p))) n converges to 
, as n tends to infinity. As a consequence, it is enough to prove that 



, [n*«J 

- ^ Trace(A 2 $ fc (p)) 



fc=[ntj_i]+l 
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converges to a strictly positive constant. A straightforward computation 
gives 

, [nU] 

lim - V Trace(A 2 $ fc (p)) 

k=[nti-i]+l 

= afaoo + a\ (1 - a^) + |a 3 | 2 + (ai + a 2 ){a 3 (3 00 + a 3 /3oo) 



(*i — * 



i-lj 



v\ + v\ + Z/| 



^l(l-^l) + ^(l-^) + ^(l-^3 2 ) 
— 21/1I/2U1U2 - 2lS 1 lS 3 ,V 1 V 3 - 2z/ 2 Z/ 3 -y 2 -y 3 • 

This means that, for every z/ 1? u 2 , v z G K, for any p > 1, the expectation 



^(n) _ X H ) + ^n) _ + ^ (z (n) _ 



converges to 



E 



PI 



where (b[ 1 \ Bf\ B^) is a 3- dimensional Brownian motion with the an- 
nounced covariance matrix. 
The polynomial 

(^f - 4- } J + - + - ^)) P 

can be developed as the sum 

S1S2 . . . 5^ 



0<Pl+P2<P 



^2 ^3 



where the summation in the last sum runs over all partitions V = {A, B,C} 
of {!,...,£>} such that |A| = pi, |S| = p 2 , |C| = P — Pi — P2, with the 



convention: 



(n) 



X. 



(n) 



(n) 



tt-i 

(n) 



*f — 1 



if j G A 

if jeB 
if jgC. 

The expectation under w of the above expression converges to the corre- 
sponding expression involving the expectation (E[-]) of the Brownian mo- 
tion (S t (1) ,S t (2) ,S t (3) ). As this holds for any z/i,z/ 2 ,z/ 3 G R, we deduce that 
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w E-p S1S2 ■ ■ ■ Sp] converges to the corresponding expectation for the Brow- 
nian motion. 

We can conclude the proof by noticing that A, can be written, modulo 
multiplication by a constant, as J2r S1S2 ■ ■ ■ S p for some p. A 

Let us discuss the class of polynomials for which Theorem 3.1 holds. In 
the particular case when the map $ is the identity map and p = 1/21 (in 
that case v i = for i = 1, 2, 3 and C = I), Biane pQ proved the convergence 
of the expectations in Theorem 13. II for any polynomial in 3m non-commuting 
variables. It is a natural question to ask whether our result holds for any 
polynomial P instead of P, or at least for a larger class. 

Let us give an example of a polynomial for which the convergence in 
our setting does not hold. Take P(X, Y) = XY. From Theorem 13.11 the 
expectation under the state uj of 

jWyW _|_ y^X^ 

converges asn-^ +00 to 2 E[S^ B\ f ]. 

Since we have the following commutation relations 

[(a x - vi I), ((T y - v 2 I)} = 2ia z , \(a y - v 2 I), (a z - v 3 I)} = 2ia x 

and 

[(a z - v 3 I), (a x - t>i /)] = 2ia y , 

we deduce that 

[Xi n \ Y t {n) ] = 2m- 1 ' 2 Z\ n) + 2itv 3 I, [Y t (n \ Z t (n) ] = 2m~ 1/2 X t (n) + 2%tv x I 
and 

\zf\ X t (n) ] = 2m- 1 ' 2 Y t (n) + 2itv 2 I. (4) 
Then the expectation under the state u of 

1 



X^Y t (n) = -[P(X,Y) + [XrX 
converges to ~E[B^ B^] + itv^ 7^ "E\B^ Bf^], if t> 3 is non zero. 

Furthermore, by considering the polynomial P(X, Y) = XY 3 + Y 3 X, it 
is possible to show that the convergence in Theorem 13.11 can not be enlarged 
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to the class of symmetric polynomials. Straightforward computations gives 
that P(X, Y) can be rewritten as 

XY^+YX~ 3 +-[X, Y}(Y 2 -X 2 )+^(Y [X, Y]Y-X[X, Y}X) + ^(Y 2 -X 2 )[X, Y] 

so the expectation w[P(Xf: n \Yt)] converges as n tends to +00 to 

E[P(S t (1) ,S t (2) )]+3ii;3t(T;?-^) 

which is not equal to E[P(S^ ,B^ )] if t> 3 7^ and \v±\ 7^ \v^\. 

In the following corollary we give a condition under which the convergence 
in Theorem 13.11 holds for any polynomial in 3m non-commuting variables. 

Corollaire 3.1. In the case when is equal to ~I, the convergence holds 
for any polynomial P in 3m non-commuting variables, i.e. for every t\ < 
£2 < • • • < t m , the following convergence holds: 



lim w 

= E 



p/ yW v( n ) y( n ) 

PCrW d(2) R (3) 



j^-(n) y(n) ^(n)\ 



R (l) R (2) R (3)> 



where {B^\ b\ 2 \ B^) 
with covariance matrix th. 



)t>o 



is a three-dimensional centered Brownian motion 



Proof: 

We consider the polynomials of the form S = J2v ^1^2 ■ ■ ■ S Pl + P2 + P3 where 
the summation is done over all partitions V = {A, B, C} of the set {1, ... , p\+ 
Vi +P3} sucn that I A \ — px,\B\ = p 2 ,\C\ — ps, with the convention: 



s 3 = {yi 



(n) 



(n) 



Y, 



(«) 



,(n) 



r(n) 



if j e A 

if J6B 

if jeC 



and N is the number of terms in the sum. 

From Theorem 13.11 the expectation under the state w of S converges to 



Using the commutation relations (jlj) with the t^'s being all equal to zero, 
monomials of S differ of each other by n^ 1 ^ 2 times a polynomial of total 
degree less or equal to (pi +p 2 +P3) — 1. It is easy to conclude by induction. 

A 
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4 Examples 



4.1 King-Ruskai-Szarek- Werner's representation 

The set of 2 x 2 self- adjoint matrices forms a four dimensional real vec- 
tor subspace of M 2 (C). A convenient basis of this space is given by B = 
{I, o~ x , o~ y , <j z }. Each state p on M 2 (C) can then be written as 

_ 1 / 1 + z x — i y \ 
^ 2 \ x + i y 1 — z J 

where x, y, z are reals such that x 2 + y 2 + z 2 < 1. Equivalently, in the basis 
B, 

p = -(I + x a x + y a y + z a z ) 

with x, y, z defined above. Thus, the set of density matrices can be identified 
with the unit ball in K 3 . The pure states, that is, the ones for which x 2 + 
y 2 + z 2 = 1, constitutes the Bloch sphere. 

The noise coming from interactions between the qubit states and the 
environment is represented by the action of a completely positive and trace- 
preserving map $ : M 2 (C) — > M 2 (C). Kraus and Choi [21 El E] gave an 
abstract representation of these particular maps in terms of Kraus operators: 
There exists at most four matrices Li such that for any density matrix p, 

$( P ) = L iP L i 

l<i<4 

with LiL* = I. The matrices Lj are usually called the Kraus operators of 
$. This representation is unique up to a unitary transformation. Recently, 
King, Ruskai et al [TUl [6] obtained a precise characterization of completely 
positive and trace-preserving maps from M 2 (C) as follows. The map $ : 
M 2 (C) — > M 2 (C) being linear and preserving the trace, it can be represented 
as a unique 4 x 4-matrix in the basis B given by 

(ii) 

with = (0, 0, 0), t e R 3 and T a real 3 x 3-matrix. King, Ruskai et al [TOl 16] 
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proved that via changes of basis, this matrix can be reduced to 



T 



( 1 \ 

ti Ai 

t 2 A 2 

\ t 3 A 3 / 



(5) 



Necessary and sufficient conditions under which the map $ with reduced 
matrix T for which |i 3 | + | A3 ) < 1 is completely positive are (see [B]) 



(Ai + A 2 ) 2 < (1 + A3) 2 - - + ^) 
(A 1 -A 2 ) 2 <(l-A 3 ) 2 -^-(t? + ^) 



1 - A 3 ± t 3 

1 - A 3 ± h 
1 + A 3 ± t 3 

2M 2 



<(l-A 3 ) 2 -i 2 (7) 



[1 - + \l + \j) - (ti + tj + u 

> 4 [\l(tj + A 2 ) + A 2 (t 2 + A 2 ) + A 2 (t 2 + A 2 ) - 2AA2A3] . (8) 



We now apply Theorem [3J] in this setting. Let $ be a completely positive and 
trace preserving map with matrix T given in (J5j), with coefficients ti, Xi,i = 
1, 2, 3 satisfying conditions §6§, ([7j) and ([8]). Moreover, we assume that |Aj| < 
1, i — 1, 2, 3. For every n > 0, 



1 + 0n(3) n (l) " 1 0n(2) 

0„(l)+i n (2) l-0n(3) 



where the sequences (0 n («))„>o, i = 1,2,3 satisfy the induction relations: 

4>n{i) = \4>n-i{i) + U. 

with initial conditions 0o(l) = %, 0o(2) = and 0o(3) = 2. Explicit formulae 
can easily be obtained. We get, for every n > 0, 



0n(l) 



^(2) = (,- I ^-)A^ I 



,(3) 



1-Ax 

*3 



A™ + 



1-Ai 

£2 



A, 



I-A3 
11 



A? + 



1-A a 



Hence, for any state p, for any n > 1, 

$ n (p) = p K 



where |AL aa . = max IAJ and 

i=l,2,3 



oo 
oo 



o(\M n max ) 



Poo 

1 — a r 



with a r 



applies with Vi 



U 



and Pa 



1,2,3. 



1-Ai 



— i 



1-A, 



Theorem 13.11 



We now give some examples of well-known quantum channels. For each 
of them we give their Kraus operators, their corresponding matrix T in 
the King-Ruskai-Szarek- Werner's representation, as well as the vector v = 
(v\, t>2, t>3) and the covariance matrix C obtained in Theorem 13. 11 It is worth 
noticing that if $ is a unitalmap, i.e. such that $(/) = /, then the covariance 
matrix C is equal to the identity matrix I3. 

1. The depolarizing channel: 

Kraus operators: for some < p < 1, 



l-pI,L 2 



IP T 



' y T 



Ip 

-<T~. 

3 



King-Ruskai-Szarek- Werner's representation: 



T 



( 1 

1-f 



\ 






1 _ i£ 
1 3 











1-f/ 



The vector v is the null vector and the covariance matrix C in this case 
is given by the identity matrix J 3 . 



2. Phase- damping channel: 



12 



Kraus operators: for some < p < 1, 



L! = ^l-pl, L 2 = y/p^ J U^,L 3 = VP^o J) 



King- Ruskai-Szarek- Werner's representation: 



T 



The vector t> is the null vector and the covariance matrix C in this case 
is given by I 3 . 



f 


1 








°\ 







1 — p 
















1 — p 





\ 











1 / 



3. Amplitude- damping channel: 

Kraus operators: for some < p < 1, 



1 



vT^p 



VP 




King- Ruskai-Szarek- Werner's representation: 

/ 1 N 

v/T^ 

VT^p 

\ t 1 -p J 



The vector v is equal to (0, 0, 1). The covariance matrix in this case is 
given by 

'10 
C = 10 10 




4. Trigonometric parameterization: 

Consider the particular Kraus operators 



Li 



v u 
cos(-)cos(-) 



/ + 



sm(-)sm(-) 
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and 



Lo = 



sm(-)cos(-) 



Ox ~ I 



cos(-)sm(-) J 



<Jy. 



King- Ruskai-Szarek- Werner's representation: 



T 



/ 1 \ 

cosm 

cost" 

y sin u sin v cos u cos v ) 



mi i / sin u sin v . rr ^ 1 
The vector v is equal to (0,0,- ). The covanance matrix 



in this case is given by 



C 



1 — COS U COS V ' 

1 
1 
1- 



with t> 3 



sin u sin v 



1 — COS U COS V 



4.2 CP map associated to a Markov chain 

With every Markov chain with two states and transition matrix given by 

p 1 — p 



is associated a completely positive and trace preserving map, denoted by 
with the Kraus operators: 







and 











y/q V 1 - Q 
Let p be the density matrix 



Vi -q, T x . \ 
— - — (/ - a z ) + —{a x - toy). 



1/1 + 2; x — i y 
2 \ x + i y 1 — z 
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where x, y, z are reals such that x 2 + y 2 + z 2 < 1. The map $ transforms the 
density matrix p into a new one given by 

= L\pL x + L;pL 2 . 

By induction, for every n > 0, 



r„ 1 - Pn 



where the sequences (p n )n>o, and (r n ) n > satisfy the recurrence relations: for 
every n > 1, 

Pn =Pn-i(p-g) + g 
and 

?"n = V g(i - g) +Pn-i(v / p( 1 -p) - Vii 1 - <?)) 

with the initial condition p = (1 + z)/2. Assumption (A) is then clearly 
satisfied with 

• q (3 



Pc 



l + q-p \ (3 1 -p 



where f3 = 
if P is a po 



q\Jp{l — p) + (1 — p)a/<?(1 — q) ■ Then, applying Theorem 13.11 
ynomial of 3m non-commuting variables, for every < t\ < t 2 < 



. . . < t m , the following convergence holds 



lim w 

n^+oo 

= E 



where (S t (1) , S t (2) , flf } ) i>o is a three-dimensional centered Brownian motion 
with Covariance matrix Ct where 

C = 

with 



1 — v\ —V1V2 

1 
— u if 2 \ — v\ 



l + q-p 



[qVp(i-p) + 0--p) ViCt-v)} 



and 



^2 



p + q—1 
l + q- p' 
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5 Large deviation principle 



Let T be a Polish space endowed with the Borel cr-field 13(F). A good rate 
function is a lower semi-continuous function A* : r— >[0, oo] with compact 
level sets {x;A*(x) < a}, a G [0,oo[. Let v = (v n ) n j oo be an increasing 
sequence of positive reals. A sequence of random variables (Y n ) n with values 
in T defined on a probability space (Q, P) is said to satisfy a Large Devi- 
ation Principle (LDP) with speed v = (v n ) n and the good rate function A* 
if for every Borel set B G 13(F), 

- inf A*(x) < liminf — logP(F„ G B) 

x&B° n v n 

< lim sup — log P(Y n G B) < - inf A* (x) . 

n V n xeB 

For every k > 1, we define 

Xfc = J(g)...<g)i'<g)cr x <g)i'<g)... 

y k = I®...®I®o-y®I®... 

z k = I <S> ■ ■ ■<%> I <S> cr z ® I <g> . . . 

where each a, appears on the k th place. 
For every n > 1, we consider the processes 

n n n 

X n = Xk, Y n — yk, Z n = Zk 
k=l k=l k=l 

with initial conditions 

A^o = Yq = Zq = 0. 

To each vector z/ = (v-±, u 2 , 1/3) G M 3 , we associate the euclidean norm \ \ = 
\Jv\ + v\ + i/| and ( ., . ) the corresponding inner product. 

Theorem 5.1. Let & be a completely positive and trace-preserving map for 
which there exists a state 



Poo = 



'00 



Poo 1 - «c 



such that for any given state p, 

$»=Poo + o(l). 
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For every v = (u ± , u 2 , u 3 ) e M 3 '* , the sequence 

( V\X n + v 2 Y n + v 3 Z n \ 

V n J n>l 

satisfies a LDP with speed n and the good rate function 



I(x) = { 



(* + h) 1o « ( 



|H±g 
IMI+<".«> 



+00 otherwise. 

where v 1 = 2 Re(/3 00 ), w 2 = -2 Im^), w 3 = 2^^ - I. 

Proof: 
The matrix 



B := l/i<7 x + Z/ 2 0"?, + v 3 a z 

has two distinct eigenvalues ±|M|. 
For every n > 0, we can write 



v 3 v x - %V 2 

vi+i v 2 —v 3 



«oo + 0n(l) A» + 0n(2) 

/3 oo + 0„(3) l-aoo + n (4) 



where the four sequences (0 n (i)) n >o satisfy </> n (i) = o(l). 

For any k > 1, the expectation of 5 in the state $ fc (p) is equal to 

Traced $ fe (p)) = ( u, v ) + e k , 

with e n = o(l). As a consequence, the distribution of B is 



l+i-r((i/,u)+e fc ) 



i-p*(-IMI). 



Using the fact that v\X n + v 2 Y n + z/3-Z n is the sum of n commuting matrices, 
we get that 

- logu? (expt(z/iX„ + v 2 Y n + v 3 Z n )) 



n 



1 n 

-$>g (elHI^dl^D + e-IHI^i-p.d^ll))) 

fc=i 
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Since e n = o(l), we obtain that 

lim - log w (exp t{v x X n + v 2 Y n + ^Z n )) 
= log ( cosh (\\v\\t) + -^jp sinh (| \v\ 

= log (cosh (| H|t)) + log (l + ^^tanh(||z/||t)). 

We denote by A(t) this function of t. 

For every i el, the function A is finite and different iable on R, then, by 
Gartner-Ellis Theorem (see [1]), the LDP holds with the good rate function 

I{x) = sup{tx — A(i)}. 

teR 

A simple computation leads to the rate function given in the theorem. 
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